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We calculate the conductance through a single quantum dot coupled 
to metallic leads, modeled by the spin 1/2 Anderson model. We adopt the 
finite-f/ extension of the noncrossing approximation method. Our results 
are in good agreement with exact numerical renormalization group results 
both in the high temperature and in the Kondo (low temperature) regime. 
Thanks to this approach, we were able to fit fairly well recently reported 
measurements £Q in a quantum dot device. We show that, contrarily to 
what previously suggested, the conductance of this particular device can 
be understood within the spin 1/2 Anderson model, in which the effects of 
the multilevel structure of the dot are neglected. 

PACS numbers: 73.23.-b, 72.15.Qm, 73.63.Kv 

The recent observation [21 El of the Kondo effect f.4_ in quantum dot 
(QD) devices has opened the possibility to control the Kondo effect experi- 
mentally, and a new exciting field of research It was shown that - below 
a Kondo temperature Tk of the order of 100 mK - the linear response con- 
ductance of a QD device, G, approaches 2e 2 jh (unitary limit) if the number 
n of electrons confined in the QD is odd, and it becomes very small if n is 
even [3]. This even/odd effect seems to be a general feature of most [H] QD 
devices. In the present paper we focus our attention on those devices for 
which Kondo anomalies appear only for odd n ^121131, and in particular 
on the experimental data reported very recently in Ref. pQ. 

In order to describe a QD coupled to its leads, we adopt the spin 1/2 
Anderson model 

H = e k4a c kfr+ £ oX! d ° d °- + Un d < [ n d i+ Y { v ^ c Lr d <^ + hx -) ■ ( 1 ) 
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Here c\ a (ck<r) creates (destroys) a conduction electron with momentum 
k, energy and spin a in one of the two leads, which we label with S 
(source) and D (drain); d\ (d a ) creates (destroys) an electron with spin 
a on the QD; is the hybridization between QD and conduction states 
(which we suppose to be k independent, \Vk a \ = Vs[d) for k G 5 or D), eq 
is the energy of a single electron localized on the QD and U is the Coulomb 
interaction among electrons on the same orbital level. In the present model 
the energy eq is not fixed, but on the contrary it is tuned by a gate voltage, 
V g , coupled to the QD through a capacitor. As a first approximation we 
assume a linear relation, Eo = aV g + const. The gate voltage (and thus Eq) 
controls the number of electrons confined in the QD at low temperature. In 
the linear response regime (Vsd "C V g , where Vsd is the source-drain bias) 
the conductance, G, may be written in a Landauer-like form jS] 

G < r - = x C ' r B im < G *( E + *>») (-1) d - < 2 > 

Here we have assumed for simplicity that the couplings to the leads are 
symmetric. The actual QD-leads coupling strength is thus T = ttN(ef)V 2 , 
where V = (Vf + V^) 1//2 , and N(ef) is the density of states per spin at 
the Fermi level |7j. In addition, / is the Fermi distribution function and 
G r (e + irj) is the retarded local Green function, i.e. the Fourier transform 
of the time-dependent function G R {t) = — i9(t)(\{d(t), eft (0)}|). In order to 
calculate the local density of states (DOS), that is p = —Im{G R (e + irj)}/Tr, 
we adopt the noncrossing approximation method in its finite-C/ extension 
(UNCA) j^j. This method allows us to explore the empty orbital, (n ~ 0, 
—Eq < — r), the Kondo (n ~ 1, T < —eo < U — T) and the doubly occupied 
(n ~ 2, — eo > U + T) regimes. Further details on the method used to 
calculate G R (e + irj) can be found in Refs. [HI HUj. 

The main purpose of the present work is to understand the experimental 
data of Ref. pQ. The results of our (new) calculations are shown in Fig. 
^b> for two different choices of parameters. The parameters of Fig. la 
are our best fit of the experimental data. We notice that the conductance 
as a function of Eq/U is directly related to experimental data because of 
the linear relation between eq and V g . The results in Fig. ^,b are in 
good agreement with previously reported exact numerical renormalization 
group (NRG) calculations (see Fig. 2 of Ref. [II]). The figure shows 
that, at high temperature (T » Tk, see e.g. T = 900 mK), there are two 
Coulomb blockade peaks at energies Eq/U ~ and Eq/U ~ —1, as observed 
experimentally by different groups ^HHII]. Each peak corresponds to the 
addition of one electron to those confined in the QD. When the temperature 
is lowered, the peaks approach each other, and the conductance gradually 
increases for n ~ 1, while it decreases for n = and n = 2. At very low 
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Fig. 1. Linear-response conductance as a function of £q/U and for different tem- 
peratures 52- ( a ) Parameters: T — 0.262 meV and U = 0.8 meV. Inset: G/Go as 
a function of T/Tk and for different choices of Eo/U. Here Go is the value of G 
at the lowest temperature that can be reached with the UNCA; Tk is such that 
G(T K ) =G /2. (b) Parameters: T = 0.120 meV and U = 0.7 meV. 

temperature (T <C Tk, see e.g. Fig. la, T = 32 mK) the Coulomb blockade 
peaks merge into a plateau at G~ 2e 2 //i, which is located in the parameter 
region for which n ~ 1. 

In the Kondo regime, the conductance is expected to be a universal 
function of T/Tk- In the inset of Fig. ^ we show the calculated universal 
curve G/Go for three different choices oieo/U in the Kondo regime (that is 
£q/U between —0.4 and —0.6). The ratio G/Go is proportional to \n(T /Tk) 
for T ~ Tk- At very low temperature (T <C T K ) G/G oc (T/Tk) 2 , as 
expected in the Fermi-liquid regime jlj and as found experimentally 

The conductance of the spin 1/2 Anderson model was calculated with 
the NRG technique in Ref. [J, and directly compared with the reported 
experimental data. Although very good agreement was reached at low tem- 
perature, it was shown that at high temperature the theoretical conductance 
strongly underestimates the experimental value. For the same choice of pa- 
rameters of the Anderson model, we find the same discrepancy (see Fig. lb 
and Fig. 2 of Ref. PQ). It was then suggested in Ref. that multilevel 
effects could be responsible of this discrepancy. Here we explored another 
possibility. We suggest that the experimental data can be fitted within the 
spin 1/2 Anderson model, provided that an appropriate choice of parame- 
ters is made. We found the best agreement for T = 0.262 meV (about two 
times larger than the one used in Fig. lb) and U = 0.8 meV. The results 
obtained with these optimal parameters are shown in Fig. la, and they are 
in excellent agreement with experiments both at high and low temperature. 
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In conclusion, we have calculated the conductance of a system made of 
a QD coupled to two leads, described by the spin 1/2 Anderson model. We 
adopted the UNCA, and we found good agreement with exact NRG results 
jllj and experiments [21 El- We have shown that recent experiments [Q, 
contrarily to what previously proposed, can be understood in the framework 
of the spin 1/2 Anderson model, in which the effects of the multi- level 
structure of the dot is neglected. Although the spin 1/2 Anderson model 
can be solved exactly by using the NRG, we believe that the UNCA is more 
suitable than NRG for extensions to realistic systems, and it could become 
an important tool to interpret experiments in which the realistic electronic 
structure of the QDs plays a crucial role (e.g. the Kondo effect in quantum 
dots for integer spin [§J). 
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